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In this article I review model-independent procedures for extracting properties of Weakly 
Interacting Massive Particles (WIMPs) from direct Dark Matter detection experiments. 
Neither prior knowledge about the velocity distribution function of halo Dark Matter 
particles nor about their mass or cross sections on target nucleus is needed. The unique 
required information is measured recoil energies from experiments with different detector 
materials. 
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1. Introduction 

Different astronomical observations and measurements indicate that more than 80% 
of all matter in our Universe are "dark" and this Dark Matter interacts at most 
very weakly with ordinary matter. Weakly Interacting Massive Particles (WIMPs) x 
arising in several extensions of the Standard Model of particle physics with masses 
roughly between 10 GeV and a few TeV are one of the leading candidates for 
Dark Matter 1 ' 2,3,4 . Currently, the most promising method to detect different WIMP 
candidates is the direct detection of the recoil energy deposited by elastic scattering 
of ambient WIMPs off the target nuclei 5,6 . The differential event rate for elastic 
WIMP-nucleus scattering is given by 1 : 



dR I poa 



dQ \ 2m x m 2 N 



F 2 (Q) 



/i(«) 



dv. (1) 



Here R is the direct detection event rate, i.e., the number of events per unit time and 
unit mass of detector material, Q is the energy deposited in the detector, po is the 
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WIMP density near the Earth, a is the total cross section ignoring the form factor 
suppression, F{Q) is the elastic nuclear form factor, fi(v) is the one-dimensional 
velocity distribution function of the WIMPs impinging on the detector, v is the 
absolute value of the WIMP velocity in the laboratory frame. The reduced mass 
m r; N is defined by 

to v ton 

™r,N = > 2 

m x + to n 

where m x is the WIMP mass and ton that of the target nucleus. Finally, 
Vmin — ay/Q is the minimal incoming velocity of incident WIMPs that can deposit 
the energy Q in the detector with the transformation constant 



(3) 



2™ r 2 N 



and u m ax is the maximal WIMP velocity in the Earth's reference frame, which is 
related to the escape velocity from our Galaxy at the position of the Solar system, 
Vcsc ^ 600 km/s. 

The total WIMP-nucleus cross section a in Eq. (1) depends on the nature 
of WIMP couplings on nucleons. Through e.g., squark and Higgs exchanges with 
quarks, WIMPs could have a "scalar" interaction with nuclei. The total cross section 
for the spin-independent (SI) scalar interaction can be expressed as 1 ' 2 

'4 



-SI 



<T,N 



Zf p + (A - Z)f n 2 . (4) 



Here m rj N is the reduced mass defined in Eq. (2), Z is the atomic number of the 
target nucleus, i.e., the number of protons, A is the atomic mass number, A — Z 
is then the number of neutrons, /( PiI1 ) are the effective scalar couplings of WIMPs 
on protons p and on neutrons n, respectively. Here we have to sum over the cou- 
plings on each nucleon before squaring because the wavelength associated with the 
momentum transfer is comparable to or larger than the size of the nucleus, the 
so-called "coherence effect" . 

In addition, for the lightest supersymmetric neutralino, and for all WIMPs which 
interact primarily through Higgs exchange, the scalar couplings are approximately 
the same on protons and on neutrons: / n ~ / p . Thus the "pointlike" cross section 
0-q 1 in Eq. (4) can be written as 

•f-(i) ^W = a'(m)" 4 , (5 ) 

where m r , p is the reduced mass of the WIMP mass m x and the proton mass to p , 
and 

: ' '"r„ .a. 2 (6) 



is the SI WIMP-nucleon cross section. Here the tiny mass difference between a 
proton and a neutron has been neglected. 
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Table 1. List of the relevant spin values of the most used spin-sensitive nuclei. More details can be 
found in e.g., Refs. 1, 7, 8, 9. 



Isotope 


Z 


J 


(S P ) 


(Sn) 


-(S P )/(Sn) 


(Sn)/(S p ) 


Natural abundance (%) 


iy p 


9 


1/2 


0.441 


-0.109 


4.05 


-0.25 


100 


' ri Na 


11 


3/2 


0.248 


0.020 


-12.40 


0.08 


100 


»»C1 


17 


3/2 


-0.059 


-0.011 


-5.36 


0.19 


76 


37 CI 


17 


3/2 


-0.058 


0.050 


1.16 


-0.86 


24 


«Ge 


32 


9/2 


0.030 


0.378 


-0.08 


12.6 


7.8 / 86 (HDMS) iu 


m T 


53 


5/2 


0.309 


0.075 


-4.12 


0.24 


100 


i^Xe 


54 


1/2 


0.028 


0.359 


-0.08 


12.8 


26 


i;ii Xc 


54 


3/2 


-0.009 


-0.227 


-0.04 


25.2 


21 



On the other hand, through e.g., squark and Z boson exchanges with quarks, 
WIMPs could also couple to the spin of target nuclei, an "axial-vector" interaction. 
The spin-dependent (SD) WIMP-nucleus cross section can be expressed as 1 ' 2 : 



,T SD 



G%m 



2 

r.N 



ill 
j 



(S p )a p + (S n )a n 



(7) 



Here Gf is the Fermi constant, J is the total spin of the target nucleus, (Sr Pta \) 
are the expectation values of the proton and neutron group spins, and a( p , n ) are 
the effective SD WIMP couplings on protons and on neutrons. Some relevant spin 
values of the most used spin-sensitive nuclei are given in Table 1. 

For the SD WIMP-nucleus interaction, it is usually assumed that only unpaired 
nucleons contribute significantly to the total cross section, as the spins of the nucle- 
ons in a nucleus are systematically anti-aligned a . Under this "odd-group" assump- 
tion, the SD WIMP-nucleus cross section can be reduced to 



,T SD 

°0 



I) 



l r,N 



111 
J 



l (p,n)\ 



(8) 



And the SD WIMP cross section on protons or on neutrons can be given as 

(9) 



a xf P ,n) = ( V ) G F W r 2 (p,n)l«(p,n)| 



V 7r 

Due to the coherence effect with the entire nucleus shown in Eq. (5), the cross 
section for scalar interaction scales approximately as the square of the atomic mass 
of the target nucleus. Hence, in most supersymmetric models, the SI cross section 
for nuclei with A > 30 dominates over the SD one 1,2 . 



2. Reconstructing the one dimensional velocity distribution function 
of halo WIMPs 

As the first step of the development of these model-independent data analysis proce- 
dures, starting with a time-averaged recoil spectrum dR/dQ and assuming that no 



a However, more detailed nuclear spin structure calculations show that the even group of nucleons 
has sometimes also a non-negligible spin (see Table 1 and e.g., data given in Refs. 1, 7, 8, 9). 
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directional information exists, the normalized one-dimensional velocity distribution 
function of incident WIMPs, fi(v), is solved from Eq. (1) directly as 11 



h{v)=M\-2Q 



dQ 



1 



F 2 (Q) 



where the normalization constant Af is given by 



Af= 2 - 

a 



dR 
dQ 

dR 
dQ 



Q =v 2/ a 2 



dQ 



(10) 



(11) 



F 2 (Q) 

Note that the WIMP velocity distribution reconstructed by Eq. (10) is independent 
of the local WIMP density po as well as of the WIMP-nucleus cross section a . 

However, in order to use the expressions (10) and (11) for reconstructing /i(i>), 
one needs a functional form for the recoil spectrum dR/dQ. In practice this requires 
usually a fit to experimental data and data fitting will re-introduce some model 
dependence and make the error analysis more complicated. Hence, expressions that 
allow to reconstruct fi(v) directly from experimental data (i.e., measured recoil 
energies) have been developed 11 . Considering experimental data described by 

Qn - h -f < Q n ,i <Qn+ b -f, i = 1, 2, ■ ■ ■ , N n , n = 1, 2, ■ ■■ , B. (12) 



Here the total energy range between Q m ; n and Q max has been divided into B bins 
with central points Q n and widths b n . In each bin, N n events will be recorded. Since 
the recoil spectrum dR/dQ is expected to be approximately exponential, in order to 
approximate the spectrum in a rather wider range, the following exponential ansatz 
for the measured recoil spectrum ( before normalized by the exposure £) in the nth 
Q— bin has been introduced 11 : 

'' m ^ V MH.,»), (13 ) 

expt, n \ / expt, Q~Q n 



dQ 



Here r n — N n /b n is the standard estimator for (dR/dQ) cxp t at Q = Q ni k n is the 
logarithmic slope of the recoil spectrum in the nth Q— bin, which can be computed 
numerically from the average value of the measured recoil energies in this bin: 



Q Qr, 



i=l 



Q„)= (y 



COtll 



k n b n 



v (14) 



Then the shifted point Q s , n in the ansatz (13), at which the leading systematic 
error due to the ansatz is minimal 11 , can be estimated by 



Qs.n Qn ~\~ * In 
Kn 



sinh(fc„6 n /2) 



k n b n /2 

Note that Q s , n differs from the central point of the nth bin, Q n . 

Now, substituting the ansatz (13) into Eq. (10) and then letting Q 
can obtain that 11 



Qs 



/l.roc (^V s ,n — OL\J Qs,n^ 



Af 



F 2 (Q s ,n) 



dQ 



lnF 2 (Q) 



(15) 



we 



(16) 



Q=Qs 
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Here the normalization constant TV given in Eq. (11) can be estimated directly from 
the data: 

-i -l 



N= 2 - 

a 



i 



VQa~F 2 {Q a ) 
where the sum runs over all events in the sample. 



(17) 



3. Determining the WIMP mass and the SI WIMP nucleon coupling 

By using expressions (10) and (11) for reconstructing the WIMP velocity distribu- 
tion function, not only the overall normalization constant Af given in Eq. (11), but 
also the shape of the velocity distribution, through the transformation Q = v 2 /a 2 
in Eq. (10), depends on the WIMP mass m x involved in the coefficient a. It is thus 
crucial to develop a method for determining the WIMP mass model independently. 

From Eq. (10) and using the exponential ansatz in Eq. (13), the moments of the 
normalized one-dimensional WIMP velocity distribution function can be estimated 
by 12 

fV (Qmax) 

(v n ) = / v n h{v)dv 

-A,(Q mi n) 

" 2 Qmit 1)/2 KQmin)/F 2 (Q min ) + (n + l)/„(Q min , Qmax) 



= a 



(18) 



1 /2 

min^ '(^3 'min) / ^ (Qmin) ^~ ^o(^?min; Qmax) 

Here v(Q) = a^/Q, G( mm . ma x) are the experimental minimal and maximal cut-off 
energies, 

r(Q min )=(^) =rie fe 1 (Q».n-Q s , 1 ) (19) 

is an estimated value of the measured recoil spectrum (dR/dQ) cxp t (before the nor- 
malization by the exposure £) at Q = Q mm , and I n (Qmin, Qmax) can be estimated 
through the sum: 

Qa n ~ 1)/2 

IniQmim Qmax) — ^ ' F 2 (Q ) ' (^) 

where the sum runs again over all events in the data set. Note that by using Eq. (18) 
(v n ) can be determined independently of the local WIMP density po, of the WIMP 
nucleus cross section cr , as well as of the velocity distribution function of incident 
WIMPs, /i(w). 

By requiring that the values of a given moment of fi(v) estimated by Eq. (18) 
from two experiments with different target nuclei, X and Y, agree, m x appearing 
in the prefactor a n on the right-hand side of Eq. (18) can be solved as 13 : 

_ v , mxiy-mx(Kii,x/V) 



rn 



X\( v n) 



Hn,xlH n .Y - y/mx/mv 



(21) 
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where 



min.X 



) + {n+ 1) J n> 



x 



^Lxrx(Q mi n,x)/F x (Q, 



.X 



1/n 



(22) 



and 1Z nt Y can be defined analogously 13 . Here n ^ 0, m^x.Y) an d F( X ,y)(Q) 
arc the masses and the form factors of the nucleus X and Y, respectively, and 
r (x.Y)(Qmin.(x,Y)) refer to the counting rates for detectors X and Y at the respec- 
tive lowest recoil energies included in the analysis. Note that the general expression 
(21) can be used cither for spin-independent or for spin-dependent scattering, one 
only needs to choose different form factors under different assumptions. 

On the other hand, by using the theoretical prediction that the SI WIMP- 
nucleus cross section dominates, and the fact that the integral over the one- 
dimensional WIMP velocity distribution on the right-hand side of Eq. (1) is the 
minus-first moment of this distribution, which can be estimated by Eq. (18) with 
n = — 1, one can easily find that 12 



PbI/pI 2 = 



(. 



1 



4^2 \£A\ 



2Q m j n r (Qmin) 
(Qmin) 



(m x + m N ) 



(23) 



Note that the exposure of the experiment, £, appears in the denominator. Since the 
unknown factor po\.f P \ 2 on the left hand side above is identical for different targets, 
it leads to a second expression for determining m x : 12 



m 



(m x /m Y ) 5/2 m Y - m x (fca,x /1Z<t,y) 



T^a,X /fca,Y - {mx/my) 

Here m^ x ,Y) o< A( X ,y) has been assumed and 



5/2 



i 



2Qm?n,X r x(Q™n,x) 
F x (Qmin,x) 



+ Io,x 



(24) 



(25) 



Remind that the basic requirement of the expressions for determining m x given 
in Eqs. (21) and (24) is that, from two experiments with different target nuclei, the 
values of a given moment of the WIMP velocity distribution estimated by Eq. (18) 
should agree. This means that the upper cuts on fi(v) in two data sets should be 
(approximately) equal . Since w cut = aVQmax, it requires that 12 



ax 

CXy 



(26) 



Note that a defined in Eq. (3) is a function of the true WIMP mass. Thus this 
relation for matching optimal cut-off energies can be used only if m x is already 
known. One possibility to overcome this problem is to fix the cut-off energy of the 



b Hereafter, without special remark all notations defined for the target X can be defined analogously 
for the target Y and eventually for the target Z. 

c Here the threshold energies of two experiments have been assumed to be negligibly small. 
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experiment with the heavier target, determine the WIMP mass by either Eq. (21) 
or Eq. (24), and then estimate the cut-off energy for the lighter nucleus by Eq. (26) 
algorithmically 12 . 

Furthermore, by combining two or three data sets with different target nuclei 
and making an assumption for the local WIMP density p 0} we can use Eq. (23) 
to estimate the squared SI WIMP coupling on protons (nucleons), |/ P | 2 - 14 ' 15 It is 
important to note that |/ p | 2 and m x can be estimated separately and from experi- 
mental data directly with neither prior knowledge about each other nor about the 
WIMP velocity distribution. 



4. Determining ratios between different WIMP nucleon cross sections 

4.1. Determining the ratio between two SD WIMP couplings 

Assuming that the SD WIMP-nucleus interaction dominates and substituting the 
expression (7) for Oq into Eq. (1) for two target nuclei X and Y, the ratio between 
two SD WIMP-nucleon couplings can be solved analytically as 16 ' 17 ' 18 

SD 



*PS ±,n 

for n ^ 0. Here I have defined 



(Sp)x ± (S p )y (Tlj,n,X /T^.J,n,Y) 
(S n )x ± (S n )Y (R-J,n,X fR-J,n,Y) 



J,n,X 



Jx \ n a . 



X 



J x + 1 / n. 



n,X 



1/2 



(27) 



(28) 



with TZ n x and Tla.x defined in Eqs. (22) and (25). 

Note that, firstly, the expression (27) for a n /a p is independent of the WIMP mass 
m x and the ratio can thus be determined from experimental data directly without 
knowing the WIMP mass. Secondly, because the couplings in Eq. (7) are squared, 
we have two solutions for a n /a p here; if exact "theory" values for Tlj. n ,(x.Y) are 
taken, these solutions coincide for 



SD 

P/ +,n 



a, 



SD 



= < 



(S P )x 

(Sn) X 

(S p )y 

(S n )Y 



for lZj, n ,x = , 
for 1Zj, n ,Y = , 



(29) 



which depend only on the properties of target nuclei (see Table 1). Moreover, it can 
be found from Eq. (27) that one of these two solutions has a pole at the middle of 
two coincident values, which depends simply on the signs of (S'rjx and (S n )Y' since 
7tj,n,x and lZ,j, n ,Y are always positive, if both of (S n }x and (S n )Y are positive or 
negative, the "— " solution (a n /a p )?. D „ will diverge and the "+" solution (a n /a p )^ D „ 
will be the "inner" solution; in contrast, if the signs of (S n )x and (S n )y are opposite, 
the "— " solution (a n /a p )^ n will be the "inner" solution. 
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4.2. Determining the ratio between two WIMP proton cross sections 

Considering a general combination of both the SI and SD cross sections given in 
Eqs. (5) and (7), we can find that 16 ' 18 



SD 

T SI 



32 \ „2 2 

— G F m 

TT J 



J+l 



r,P I 



(S p ) + (S n )(a n /a p ) 



XP 




, (30) 



where af p given in Eq. (9) has been used and 



(S p ) + (S n )(a n /a p ) 



1 2 



C = t( J -± 

Then the expression (1) for the differential event rate should be modified to 



(31) 



(dR\ 



= £A 2 \ 



expt, Q=Q„ 



2m x m 2 iP 



F SI (Q m in) + -§?- C P F SB (Q 



f(Qmax) 



h{v) 



dv , 



(32) 



u(Qmin) 

where I have used Eq. (5) again. Now by combining two targets X and Y and assum- 
ing that the integral over the WIMP velocity distribution function in Eq. (32) esti- 
mated by Eq. (18) for each target with suitable experimental maximal and minimal 
cut-off energies should be (approximately) equal, the ratio of the SD WIMP -proton 
cross section to the SI one can be solved analytically as 16 ' 17,18 



SD 
XP 
• SI 

XP 



-^SI y(Qmin,y)(7^m,x/7^m,y) 



.X 



) — Cp,YF ST) y (Qmin, y)(^m,X /7^m,y) 



where I have assumed m 



(X,Y) 



oc A(x,y) an d defined 



n 



m,X 



rx(Qmin,x) 



(33) 



(34) 



Similarly, the ratio of the SD WIMP-ncutron cross section to the SI one can be 
given analogously as d : 

^ SL) FsiY (Qmin, y)(^m,x/^m,y) — -^fl X {Qmin,X ) 



X" _ 



c^p C ILt xF ST)X (Q m i IL! x) — C n y F^ TyY {Q ra i n ,y)(TZ mt x l^m.y) 



with the definition 




(S P )(a p /a n ) + (S n ) 



A 



(35) 



(36) 



Note here that one can use expressions (33) and (35) without a prior knowledge 
of the WIMP mass m x . Moreover, a x(p n)/ a xp are lunc ti° ns 01 omv T^-m,(x,Y), 
or, equivalently, the counting rate at the experimental minimal cut-off energies, 



d Here I assumed that a^} n ~ cr|p. 
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r (x,Y)(Qmin,(x,Y))i which can be estimated with events in the lowest available en- 
ergy ranges. 

On the other hand, for the general combination of the SI and SD WIMP- 
nucleon cross sections, by introducing a third nucleus with only the SI sensitivity: 
(S P )z = (S n )z — 0, i.e., Cp.z = 0. The a n /a p ratio can in fact be solved analytically 



as 



16,17,18. 



SI+SD 



- {Cp,XS n /p,X ~ Cp,ySn/pj) ± y/Cp,XCp,Y 



,X — b n/p,Y 



Here I have defined 



C P. xS n/p,X C P^ S n/p,y 



■ (37) 



C P,X = 77 



4 ( Jx + 1 



ix 



(S P )x 



Fsi,z(Qmin,z) 



K 



m,Z 



Fgi Y (Qmin,y) 



^SD,X 



), (38a) 



Cp.Y 



Jy + 1 
Jy 



p)Y 



Ay 



Fsi,z(Qmin,z) 



rn.X 



m,Z 



- F Sl,x(Qrain,x) 



-^SD,r(Qmin,y) ', 



(38b) 



and s n / P ,x = (S n )x/(S p )x- Note that, firstly, (a n /a p )| I+SD and c p ^ x ,Y) given in 
Eqs. (37), (38a), and (38b) are functions of only r^ x ,Y,z)(Qmin,(x.Y,z)), which can 
be estimated with events in the lowest available energy ranges. Secondly, while the 
decision of the inner solution of (a n /a p )j l P n depends on the signs of (S n )x and 
(S n )y, the decision with (a n /a p )j. I+SD depends not only on the signs of s n / P ,x = 
(S n )x/(Sp}x and s n / p _ Y = (5 n )y/(5' p )y, but also on the order of the two targets. 

Moreover, since in the expression (33) for the ratio of two WIMP-proton cross 
sections there are four sources contributing statistical uncertainties, i.e., C Pi (x,y) 
and lZ m ^x,Y), in order to reduce the statistical error, one can choose at first a 
nucleus with only the SI sensitivity as the second target: (S p )y = (<S' n )y = 0, i.e., 
Cp.y = 0. Then the expression in Eq. (33) can be reduced to 18 



SD 

xv 



^SI y( < 3min,Y)(7£m,x/7£m,r) 



-^SI,x('5n 



,x) 



C 



P.X-fSD.xlQmin,* 



:) 



(39) 



Secondly, one chooses a nucleus with (much) larger proton (or neutron) group spin 
as the first target: (S p )x ^> (S n )x — 0. Now C p .x given in Eq. (31) becomes 
(almost) independent of a n /a p : 



--p,x 




(S P ) 



p)X 



(40) 
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5. Summary and conclusions 

In this article I reviewed the data analysis procedures for extracting properties of 
WIMP-like Dark Matter particles from direct detection experiments. These meth- 
ods are model-independent in the sense that neither prior knowledge about the 
velocity distribution function of halo Dark Matter nor their mass and cross sections 
on target nucleus is needed. The unique required information is measured recoil 
energies from experiments with different target materials. 

Once two or more experiments observe a few tens recoil events (in each exper- 
iment), one could in principle already estimate the mass and the SI coupling on 
nucleons as well as ratios between different cross sections of Dark Matter parti- 
cles. All this information (combined eventually results from collider and/or indirect 
detection experiments) could then allow us to distinguish different candidates for 
(WIMP-like) Dark Matter particles proposed in different theoretical models and to 
extend our understanding on particle physics. 
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